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ABSTRACT 

We compute the phase lags between the radial velocity curves and the light curves A$i = <j^ r - (f>™ ag for 
classical Cepheid model sequences both in the linear and the nonlinear regimes. The nonlinear phase lags 
generally fall below the linear ones except for high period models where they lie above, and of course for low 
pulsation amplitudes where the two merge. The calculated phase lags show good agreement with the available 
observational data of normal amplitude Galactic Cepheids. The metallicity has but a moderate effect on the 
phase lag, while the mass-luminosity relation and the parameters of the turbulent convective model (time- 
dependent mixing length) mainly influence the modal selection and the period, which is then reflected in the 
Period - A$i diagram. We discuss the potential application of this observable as a discriminant for pulsation 
modes and as a test for ultra-low amplitudes (ULA) pulsation. 
Subject headings: Cepheids — instabilities — stars: oscillations 



1. INTRODUCTION 

The relative phase between the luminosity and the velocity 
curves in classical pulsating variables was a puzzle in the early 
days of variable star modeling. Because, overall, the pulsa- 
tions are only weakly nonadiabatic it was expected that the 
maximum brightness should occur at maximum compression, 
i.e., minimum radius, whereas in reality it is observed to occur 
close to maximum velocity. (To avoid confusion we note that 
in this paper we use astro-centric velocities, u = dR/dt.) The 
luminosity thus has a ~ 90deg phase lag compared to the one 
expected adiabatically. 

As it became possible to make accurate linear and nonlin- 
ear calculations of the whole envelope, including in particular 
the outer neutral hydrogen region, agreement betwee n mod- 
eling a nd observation was achieved. However, it was ICastoj 
(1968) who first provided a qualitative physical understand- 
ing of this phenomenon. He pointed out that during the com- 
pression phase, for example, the hydrogen partial ionization 
front moves outward with respect to the fluid, and energy is 
removed from the heat flow to ionize matter. Because of this 
temporary storage of energy - as in the charging of a capacitor 
- the H acts acts like a low-pass filter, causing a phase lag of 
close to 90 degrees. However, the exact value of the phase 
lag depends sensitively on details of the stellar model, and 
can only be obtained by detailed simulations. An understand- 
ing of the phase lag involves physics both in the linear and 
the nonlinear regimes. This characteristic makes it an ideal 
benchmark to test existing hydrocodes against observational 
constraints. 

Cepheids are known to have humps or bumps on the light 
curves, especially near maximum light. Rather than define 
the phases with respect to maximum light, or some other cho- 
sen reference point, it is therefore desirable to define them in 
terms of a Fourier sum: 



/(f) = 5Z A * sm(kajt + (/) k ) 



(1) 



The phase lag A$i ={(f)\ T — ^>™ g ) is thus defined as the dif- 
ference between the phases of the fundamental components 
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of the Fourier fits of the radial velocity (V r ) and magnitude 
curves. The two time sequences are of course reduced to the 
same time origin (epoch). 

Good Galactic Cepheid light curve data have been available 
for a long time, wh ile radial velocity data are more recent. 
Qgloz aef al. I ((2000) have published the phase lags for a set 
of fundamental (F) and first overtone (01) Cepheids and have 
found that for the F Cepheids A$i is largely independent of 
pulsational period, with an average value of -0.28, but that 
for the Ol Cepheids it decreases with period from -0.24 to 
-0.71. 

In this paper we describe the results of an extended sur- 
vey of the linear and the nonlinear (full amplitude) Cepheid 
pulsation models, and we revisit the phase lag problem. In 
our computations we make use of the state-of-the-art Florida- 
Budap est hydrocode (lYecko et al. Ill 9981: iKollath et al. 1119981 
2002) which includes a model for turbulent convection. 

Our goal is to compare the modeling results with a turbu- 
lent convective hydrocode for the phase lag with the available 
observational data. The last si milar systematic phase lag in- 
vestigation was carried out by Simon & Davisl (|1983l) using 
early radiative hydrocodes and opacities. Furthermore, the 
limited number of their models, and the quality of the obser- 
vations that were available at the time justify a new look at the 
problem. 

2. LINEAR AND NONLINEAR CEPHEID MODELING 

With the usual convenient, but excellent assumptions of a 
chemically homogeneous envelope and a rigid, nonpulsating 
core, we can then specify a Cepheid envelope by three pa- 
rameters for a given composition: effective temperature (r e t), 
luminosity (L) and stellar mass (M). The composition is usu- 
ally specified by the H content X, the metallicity Z and the He 
content Y=l-X-Z. 

Cepheids evolve on tracks that are provided by stellar evo- 
lution calculations L = L(t;M) and T e f = T e s(t;M), and the vast 
majority of Cepheids are expected to be found on the sec- 
ond, blueward penetration of the instability strip (IS). Since 
the Cepheid tracks are more or less horizontal (fixed L) it is 
customary and expedient (but not necessary) to identify these 
tracks by an M-L relation, L = L(M ; X ,Z). This reduces the 
number of stellar parameters by one, and we can thus make 
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FIG. 1 . — The Resonance, Girardi and Bono mass-luminosity relations for 
Galactic, LMC and SMC metallicity. 

one-parameter sequences of models for a given mass M and 
composition, with T e $ as the single parameter. The selected 
M-L relation however affects the whole picture from modal 
selection, to mass-period relation, to phase lags. 

The M-L relations are usually obtained from evolution- 
ary calculations. The problems with evolutionary tracks 
producing Cepheids and the difficulty of deriving mass- 
luminosity relations are well known and are addressed in 
iBuchler & Szabol (120071) . In this paper we make use of three 
different M-L relations. Our first M-L relation is a fit to the 
evolutionary tracks w ith a non zero overshooting p arameter of 
iGirardi ef at] (120001) . Second. IBotio et al. I (120001) give an M- 
L relation that is a fit to their evolutionary tracks without over- 
shooting. Third, one can also bypass evolution calculations 
and combine observational data with pulsatio nal theory. This 
has be en done for the SMC and the LMC by lBeaulieu et al. I 
d2001l) who derived stellar masses from the observed magni- 
tudes, colors and periods, but reddening errors prevented the 
extraction of a very tight M-L relation. An alternative is to 
pinpoint the P2: Pq=1 :2 resonance center from the behav- 
ior of the Fourier phase 02 1 as a function of period, and then 
use linear modeling to match the resonance condition and thus 
tie down a 'Resonance' M-L relation, which is assumed to be 
linear: LogL = a + 3.56LogM, where a is 0.7328, 0.96864, 
1.1552 for Z = 0.020, 0.010 and 0.004, respectively. We refer 
to it as the Resonance M-L relation. (In this paper we label 
the modal periods with numeral subscripts, i.e., for F, 1 for 
Ol, etc.) 

The three M-L relations are displayed in Fig. Q] For 
Z=0.020 the Resonance M-L relation runs close to the Bono 
relation in the low mass regime, while the Girardi relation 
predicts much lower luminosity. In the high mass range 
(M > 5.8M Q ) the Resonance M-L relation runs between the 
Girardi and Bono formulae for this metallicity. This charac- 
teristic has led us to choose the Resonance M L relation as the 
reference throughout this paper. The sensitivity to the form of 
the M-L relation will be discussed in Sec. 14.31 

The mass range of our survey extends from a lower M of 
2.0M© to an upper M that we vary with the metallicity and ob- 



TABLE 1 

Turbulent convection parameters 





Old 


Otc 


at s 


ot v 




Oir 


a,, 




A 


2.177 


0.4 


0.4330 


0.12 


0.001 


0.4 


0.0 


1.5 


B 


2.177 


0.4082 


0.4082 


0.25 


0.0 


0.0 


0.0 


1.5 



servational constraints, i.e., the observed longest period pulsa- 
tion, and is 9.5M , 7.5M Q and 6.5M for Galactic, LMC and 
SMC variables, respectively. The mass increment is 0.5 M Q 
between the sequences. The luminosity is obtained from the 
M-L relation. The temperature range covers the whole classi- 
cal IS, and subsequent models are computed 100 K apart. 

In our calculations we employ the Florida-Budapest hy- 
drocode (lYecko et al. I fl998l: kollafh et al.] fl998l l2002h to 
produce equilibrium Cepheid models, to perform their lin- 
ear stability analysis and to compute full amplitude pulsa- 
tions. This code incorporates a time-dependent mix ing length 
model f or turbulent convection (|Yecko et aOl998l) . The G91 
OPAL (Iglesias & Rogers 1996) with solar element abun- 
dances dGrevesse & Anders! 1 99 ll) and lAlexander & Fer guson 
d 19941) opacities are used. The equation of state and the opac- 
ity are given in terms of the usual composition variables, X, Y 
and Z. 

The values of the turbulent convective (a) parameters in the 
code affect most of the observables of Cepheid models, such 
as the width of the IS, the shape of the light and radial ve- 
locity curves, the amplitudes, the locations of the resonances, 
the maximum overtone period, just to name a few. It has been 
found that very different parameter sets can give surprisingly 
similar results. But a global calibration of the turbulent con- 
vective parameter set has not been performed yet, because of 
the huge computational task. Therefore we restrict ourselves 
to only two sets, which are presented in Table Q] The 8 pa- 
rameters are dimens ionless, of order of unity, but they obey a 
scaling invariance dKollath et al Il2002l). For a comparison to 
other mixing length models see lFeuchtinger et aT\ d2000t) . 

The linear code provides the radial displacement, veloc- 
ity, temperature and turbulent energy eigenvectors for all the 
radial pulsation modes at any point, and in particular in the 
photosphere. From these we can compute the luminosity and 
magnitude eigenvectors, and then the phase lags A$i. In this 
paper we are interested in periodic F, Ol and 02 pulsations 
(limit cycles). 

The equilibrium models are initiated with a velocity kick in 
one of the eigenmodes and the desired full amplitude (limit 
cycle) pulsations are computed with our hydrocode. The 
regime where the F and Ol modes are simultaneously linearly 
unstable is well known to split into two types, depending on 
the stellar parameters, namely a double-mode regime, where 
both modes achieve full amplitude, and a hysteretic regime - 
the 'either-or' regime in the pulsation jargon - where F and 
Ol limit cycles coexist, and the initial kick determines which 
of the two limit cycles is attained. A similar split occurs for 
the Ol and 02 modes. (In some narr ow regions the s ituatio n 
can be even more complicated, e.g., Kolla thef al. I d2002l) h 
The double-mode pulsations (also called beat pulsations) are 
not of interest to us here, but are taken up in a subsequent pa- 
per. The desired phase lags A$i are obtained from the mag- 
nitude and V r curves of the computed limit cycles through an 
8th order Fourier fit (Eq.Q]). 
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FIG. 2. — Schematic HR diagram showing the linear instability strips of 
the fundamental (F), first (Ol) and second overtone (02) modes. At the solid 
lines the linear and nonlinear edges coincide. 

3. MODAL SELECTION - THE BIFURCATIONS 

First we recall the general topography of the linear IS 
shown schematically in an HR diagram on Fig. [2] In the 
domain of astrophysical parameters that are relevant for 
Cepheids, only the F, Ol and 02 modes are unstable. As to 
the definitions of linear and nonlinear edges we refer to the 
literature, e.g., (IStellingwerj|1975HBuchler & Kovacsll 19861: 
Szabo & Buchle j[2007h and in particular Figs. 4 and 5 in 
Feuchtinger et al.\ yOOO) whose notation we will use here, 
namely L for linear and N for nonlinear, with BE for the blue 
edge and RE for the red edge. Fig.|2]shows that the region of 
linear instability of the F mode extends to very high luminosi- 
ties. That of the Ol mode resembles a sugar cone, thus with a 
maximum luminosity, which overlaps with the F IS. Similarly 
the 02 IS sugar cone, with a yet lower maximum luminosity, 
overlaps with both the 01 and F ISs. Along the solid lines the 
linear and nonlinear edge s coincide, and the pulsa tion ampli- 
tudes vanish there (See Buchler & Kollath (2002) for a more 
thorough discussion of how the behavior at the edges of the IS 
is affected by stellar evolution). The dashed parts of the linear 
instability boundaries are not real edges because an evolving 
star that is pulsating in a given mode undergoes a switch (bi- 
furcation) to another type of pulsation before it can reach the 
dashed line corresponding to the original pulsation mode. The 
substructure in the topogr aphy of the IS will be d iscussed in a 
paper on Beat Cepheids dSzabo & B uchler 20071). 

For the purposes of this paper we note that because the star 
pulsates with infinitesimal amplitude near the solid edges, the 
linear and nonlinear phase lags therefore approach each other. 
All other mode switching occurs with finite amplitude and the 
linear and nonlinear A$i differ therefore at the bifurcation. 
This behavior will clearly be visible in the subsequent figures. 
Thus for example, a star on a blueward horizontal track at 
high luminosity in Fig. [2] starts to pulsate with infinitesimal 
amplitude when it crosses the FRE (the solid line on the right) 
and stays in this mode until it reaches the FBE (solid line) on 
the left where its pulsation amplitude goes to zero. The linear 
and nonlinear A$i should thus be equal to each other both at 
the FRE and at the FBE. In contrast, a star on a blueward track 
at mid-magnitude in Fig. [2] similarly starts to pulsate when it 
crosses the FRE (the solid line on the right), but it bifurcates 
to some other pulsational state before it reaches the LFBE (the 
dotted line corresponding to the F mode). The amplitude of 
the pulsation does not go to zero at the switching which has as 
a consequence that the linear and nonlinear A$i are different. 
Finally, we note that because of the M-L relation, the ordinate 
in the schematic Fig. [2] could also represent the stellar mass 
M. 



In the two hysteretic regimes, where the pulsation is peri- 
odic, it can be in either the F or the Ol mode, or it can be 
in either the Ol or 02 mode, respectively. In this regime we 
include the corresponding model in both the F and the Ol set 
of pulsators. The double-mode pulsations of two types that 
have been encountered, namely F&Ol and 01&02, are ig- 
nored here because of their mu lti-periodic nature. The y will 
be discussed in a separate paper (Szabo & Buchlerl2007l) . The 
temperature range of the double-mode region is typically less 
then \Q0K, therefore the exclusion of these models has little 
effect on the results of this paper. 

4. LINEAR AND NONLINEAR PHASE LAGS 
4.1. Results 

The general picture of the theoretical linear and nonlin- 
ear phase lags is displayed in Fig. [3^. These results are ob- 
tained by using the Resonance M-L relation, the turbulent 
convective parameter set A (Table [TJ, and a Galactic metal- 
licity (Z=0.020). Below we investigate how each of these 
three choices affects the results. The phase lags are plot- 
ted as a function of the logarithm of the period of the ex- 
cited mode (The models in the 'either F or Ol' regime appear 
twice, once under Po and once under P\). Our models form 
one-parameter sequences at a given M, and thus L = L(M), 
with r e f as the variable parameter. The masses of the se- 
quences run from 2.5M Q to 9.5M Q , with a 0.5M Q increment, 
and are labeled by the small (integer) numbers. The F, Ol and 
02 groups start with the same masses on the left. The dotted 
lines depict the linear phase lags of the linearly unstable mod- 
els, while the thick lines correspond to the nonlinear phase 
lags, i.e., those of the limit cycles. (Some inter/extrapolation 
has been done near the edges.) 

In Fig. [3^ one can identify the three separate clusters that 
correspond to the F, 01 and 02 radial pulsation modes of 
the sequences of models. For high masses, the nonlinear se- 
quences follow the linear ones of the same mass. At a certain 
mass Ol pulsation appears, thus the nonlinear F curves stop 
before reaching the FBE, and the sequence continues in the 
corresponding Ol curve. In the 'either-or' regions, the same 
model contributes to two pulsational modes. The same thing 
happens where the sequences change from Ol to 02 pulsation 
at even lower mass. Therefore the F and Ol sequences seem 
to be increasingly incomplete toward the short periods at low 
mass. At the lowest mass regime, the models reach the blue 
edge pulsating in the 02 mode. 

As expected, the F mode pulsations cover the broadest 
range in period and in mass, as well as in phase lag. We 
remark in advance that when comparing the theoretical val- 
ues of the phase lags to observations the most relevant part of 
these arches are the broad vicinity of their maxima, because 
the long 'legs' on the high period (right) side of the sequences 
are associated with low pulsational amplitudes, near the FRE. 
We will return to this issue. Along the sequences the dots 
mark 100K intervals which allows one to see that the mod- 
els are not equidistant along the A$i in Fig. [3] but that they 
crowd toward the center of the IS and are increasingly spaced 
out toward the FRE. The FRE is the lower envelope of the F 
sequences on the figure, corresponding to infinitesimal pulsa- 
tional amplitude. The spacing is increasing toward the low 
mass models, too. The maximum F phase lags of the models 
thus lie close to zero, while the bulk of our models are in the 
-0.1 to -0.4 range. 

The phase lags for the Ol and the 02 pulsations of the se- 
quences are located lower on the diagram to the left. With 
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FIG. 3. — a.) Linear and nonlinear phase lags of sequences of equal mass Galactic Cepheid models. The mass labels appear on top. The dashed lines denote 
the linear phase lags of linearly unstable models, the thick lines denote the nonlinear (limit cycle) phase lags. The triangles along the model sequences denote 
100K intervals in r c j. The Resonance M-L relation was used along with the a parameter set A (Table[T}. b.) The difference between the nonlinear and linear 
phase lags (A(A<I > i)) for the same models, as a function of the pulsational period. The Ol and 02 models are shifted vertically by -0.2 and —0.1, respectively 
for clarity. See color version of the figures in the online edition. 

the chosen set A of a parameters the maximum Ol period is 
~ 7 d , and the maximum 02 period is <~ l d . The Cepheids pul- 
sating with the maximum Ol period are close to the LOBE, 
and therefo re low amplitudes a r e exp ected to hamper their dis- 
covery (see Feuchtinger et al. ( 200(J for a similar comment). 
The same argument applies for the maximum 02 period, al- 
though no Galactic 02 Cepheids have been found to date. 

We need to point out that this maximum Ol period is sen- 
sitive to the turbulent convective parameters that we use. This 
will be discu ssed in S ec. l4~4l and a thorough calibration will be 
attempted in Szabo & Buchler (2007). We note however, that 
our phase lag results do not seem very affected by changes in 
as as long as we use a parameter combination that reproduces 
the maximum Ol period. 

4.2. Resonances 

It is well known that internal resonances cause structure in 
the Fourier coefficie nts vs. period or some other parameter, 
(e.g jBuchlerl (119931) ), Thus the light curves of RR Lyrae who 
have no low-order internal resonances exhibit very smooth 
Fourier decomposition coefficients (<^2i,^2i> etc), whereas the 
Cepheids which span a much wider range of M, L and T e $ en- 
counter a number of such resonances durin g their evolutio n 
through the IS (see for example Fig. 4 in Buchlerl dl997l) ). 
We will see that these resonance characteristics play a role in 
this study. This may come as a surprise because one might 
have thought that resonances affect for light curves and ra- 
dial velocity curves and that their effect would balance out in 
A$i = cj)\ r — <p™ as . This is however not the case because the 
resonance induced distortions are different for light and for 
radial velocity. 

As we mentioned earlier, the differences between the non- 
linear and linear phase lags, A(A$i), vanish at the IS edges. 



This is depicted on Fig. |3^, especially for the F and Ol se- 
quences. The |A(A$i)| are very small for the 02 pulsators 
which have very low amplitudes. Generally, one expects that 
the larger the amplitude the larger the difference between the 
linear and nonlinear quantities such as the phase lag. This 
simple argument may break down where there are resonances 
which distort the light and radial velocity curves. 

The largest |A(A$i)| are found around P = I0 d for the F 
pulsators, but it is not immediately clear whether this is re- 
lated to the well-known Pi : Pq = 1 : 2 resonance, or whether 
it is just a broad minimum occurring there by coincidence, 
as part of a general trend running through the whole mass 
range. To resolve the issue, we display the A(A$i) values 
on Fig. |3J). The Ol and 02 sequences are shifted vertically 
for clarity. The A(A$i) are positive for the highest mass val- 
ues, and tend to be negative for lower masses. For a given 
mass sequence its absolute value is largest toward the middle 
of the instability strip. We point out that this figure is very 
similar for the A(A$i) for LMC and SMC metallicities that 
are discussed in Sec. 14.51 

The A(A$i) tend to be more negative in the vicinity of 
the resonance, and they do not extend to the positive region. 
On the low period side A(A$i) changes abruptly due to the 
distortion of the light and velocity curves that is caused by 
the resonance. Based on these findings, we conclude, that the 
P2 : Pq = 1 : 2 resonance has a noticeable effect on the phase 
lags. Namely, an additional dip appears on the upper envelope 
of the arched nonlinear phase lags which is very pronounced 
around the resonance center. The Log P - A(A$i) diagram 
shows, that the resonance related distortion plays an equally 
important role in shaping the phase lag curves as the ampli- 
tudes themselves. The same resonance behavior will be no- 
ticed in the LMC and the SMC metallicity models around the 
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FIG. 4. — Linear phase lags for Galactic metallicity (Z=0.020) with three M-L relations: a.) Girardi, b.) the Resonance M-L relation, and c.) Bono, for the a 
parameter set A. 



apposite resonance periods confirming our argument. 

Similarly, we believe that the surprisingly large | A$i | that 
are found for the Ol modes of the 5.5 M© sequence are 
associated with the P4 : Pi = 1 : 2 resonance near an over- 
tone period Pi = 4?2 period. This resonance causes the 'ei- 
ther or region' to broad en to a width of several hundred K 
( Szabo & Buchlerl 120071) . which is in turn reflected in the 
much longer Ol phase lag curve. The much shorter 6M Q 
curve is still affected by the resonance, as its strange shape 
and the large A(A$i) indicate. 

The effect of these resonances can also be witnessed on the 
amplitude behavior which is discussed further in Sec. [5] 

4.3. Effects of the M-L relation 

Fig.EJcompares the linear A$i for the three M-L relations 
that we have described in §2. In order to avoid overcrowding 
we plot only the linear phases lags. The overall difference 
between the nonlinear and linear lags are found to be similar 
to what we have seen in Fig.|3^ for the Resonance M-L . 

One immediately notices that the effect of the three M-L re- 
lations is most prominent in the high mass regime, where the 
highest periods are quite different. The Bono M-L relation 
has lower luminosity than the other two, and it would require 
a much higher mass (> 1OM ) to reach the 100 rf fundamental 
period limit to which the 9.5M Q sequence rises for the Reso- 
nance M-L . On the low-mass and low luminosity end the 1 ?0 
F period occurs for 3.50, 3.00 and 3.25M Cepheids, for the 
Girardi, Resonance and Bono M-L relations, respectively, as 
one would expect from Fig. Q] 

With the Girardi relation the A$i values of the models turn 
upward (especially on the low-temperature (right) side). The 
same trend is visible with the Resonance M-L relation but 
does not occur with the Bono one, again consistently with the 
behavior of the M-L relations. 

We conclude that the effect of the M-L relations on the 
A$j -logP diagram occurs indirectly through a modification 
of the mass-period relation. 

Generally the observed phase lags, which are shown in 
Fig. [4] and the mass-period distribution show good agreement 
with the models obtained by the Girardi and the intermedi- 
ate M-L relations, but disagree with the Bono relation at high 
mass. 

4.4. Effects of the a Parameters 



The two sets of convective a parameters that we use here 
are given in the Table. The main difference lies in the amount 
of turbulent viscosity (a u ) which is higher for set B. Set A 
features a small amount of overshooting (a,) and Peclet- 
correction (a r ), which accounts for the decrease of convective 
efficiency when radiative losses are important. Despite the 
smaller turbulent viscosity in set A, the models computed with 
these parameters show lower luminosity and velocity ampli- 
tude by a factor of two compared to models in set B. The 
smaller a„ is overcompensated by the remaining parameters. 

Fig. [5] compares the linear and nonlinear phase lags com- 
puted with set A ( a.) panel) and set B ( b.) panel) convective 
parameter sets. The Resonance M-L relation and a Galactic 
metallicity (Z = 0.020) have been used in both. 

The shape of the phase lag curves is similar on the two pan- 
els, but the range of the phase lags is smaller for set B. The dif- 
ference has to do with the lower pulsation amplitudes that are 
found for set B. Another noticeable feature is the smaller dif- 
ference between the linear and nonlinear phase lags for panel 
B. The underlying cause is again the too low amplitude with 
the set B parameters, because this difference depends on the 
nonlinear effects, i.e., the magnitude of the amplitude. 

The convective parameters affect the modal selection. Our 
main concern has been to reproduce the observed maximum 
Ol period (7?57 for V440 Per). We find that the longest 
period among the nonlinear overtone models is 7?7 for the 
set A models, and a very small value for the set B parame- 
ter choice (less than 3 days) that is unacceptable on observa- 
tional grounds. The difference in modal selection manifests 
itself in the 02 pulsation, too. For set A, there is (pure) sec- 
ond overtone pulsation for mass lower than 3.5 M . However 
the radial second overtone mode is linearly unstable only for 
masses below 2.5 Mq in set B, and no 02 pulsation was found 
in the models in the hydro grid in this case. This means, that 
if 02 pulsation exists there, it would occupy a rather narrow 
temperature range. 

Our experience is that, as long as they are capable of repro- 
ducing the observed mode selection patterns, the turbulent a 
parameters do not have a major effect on the phase lags. 

4.5. Effects of the Metallicity 

We have adopted here the following composition pairs (X, 
Z ) = (0.700, 0.020), (0.716, 0.010), (0.726, 0.004) as 'canon- 
ical' values when computing Galactic, LMC, and SMC mod- 
els, respectively. 
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FIG. 5. — Linear and nonlinear Galactic phase lags computed with two different a parameter combinations a.) set A and b.) set B. and using the Resonance 
M-L relation. Dashed lines: linearly uns table models, thick l i nes: nonlinear (limit cycl e) phase lags. The observed phase lags with error bars are overplotted for 
F (squares) and Ol (triangles) Cepheids (Ogloza et al. 2000; Moskalik & Ogloza 2000). See color version in the online edition. 
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FIG. 6. — The effect of metallicity on the phase lags. LMC and SMC models are shown on panel a.) and b.), respectively. Dashed lines: linearly unstable 
models, thick lines: nonlinear (limit cycle) phase lags. The results are obtained with the Resonance M-L relation and the a parameter set A. 
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In Fig. [6] we display the linear and nonlinear phase lags 
for the LMC ( a.) panel) and SMC ( b.) panel) metallic - 
ities. For short periods the phase lag is insensitive to the 
metallicity, but for high periods (i.e., high masses) lowering 
the metallicity causes the theoretical F phase lags to shift up- 
ward, and even causes a reversed sign. The same trend can 
be seen for Ol models, although to a lesser extent. We note 
that the P2 : Po = 1 : 2 and P4 : Pi = 1 : 2 resonances mentioned 
in Sec. 14.21 have very similar effects on the phase lag of the 
lower metallicity model sequences. Unfortunately no obser- 
vational extragalactic A$i are available to date on which we 
could test our predictions. 

The maximum period of the possible 02 pulsation de- 
creases with decreasing Z, but observationally, a number of 
02 and 01&02 double mode Cepheids are known in the 
Magellanic Clouds in contrast to the Galaxy. This does not 
cause a problem though, because the evolutionary tracks also 
shift with Z. Also, our maximum predicted 02 periods are 
close to maximum 02 perio d reported for Q1&Q2 double 
mode Cep heids in the LMC ( Soszyns kFef al. 1 1200 lb and in 
the SMC dUdalski et al. |[l999l) . 

So far we have only considered a solar mix where the 'met- 
als' are lumped collectively into Z. However, it has been 
found th at deviations from the stan dard solar chemical com- 
position dGrevesse & Andersll 19911) play a nonnegligib le role 
in beat Cepheid period ratios (Buchler & Szabo 2007). We 
have therefore also computed the linear phase lags for a model 
sequence with th e same opacities as used and described in 
iBuchler & Szabol (|2007), namely Fe group element number 
densities arbitrarily reduced to 75% at fixed Z. The phase lags 
are found to be largely insensitive to the chemical buildup, 
the difference is largest for high mass models, being less than 
0.01, and in general it is much less. 

After having demonstrated that the phase lags are moder- 
ately sensitive to metallicity and that they implicitly depend 
on the M-L relation and the turbulent convective parameters, 
we examine other possible physical and numerical effects on 
the Cepheid phase lags. For these tests we have computed 
only linear models, because the relationship between nonlin- 
ear and linear phase lags is not strongly affected, and this will 
therefore not significantly alter our conclusions. 

4.6. Stellar Rotation 

In this section we consider the effect of stellar rotation be- 
cause it has the effect of rearranging the stellar structure, and 
can cause slight changes in the phase lags. However, since 
we do not have the 3D code to conduct a proper calculation 
of the effects of rotation, we resort to an expedient subterfuge 
which we expect to be accurate enough as an order of magni- 
tude estimate. As in Buc hler & Sza bo (2007) and references 
therein we incorporate in our ID hydrocode a radial pseudo- 
centrifugal force of the form (F = u 2 r) which, if anything, 
is likely to overestimate the effects of rotation. We find that 
even up to 20 km s" 1 rotation, the phase lags are indistinguish- 
able from the nonrotating ones, except perhaps for the high- 
est masses (>9M Q ), where the relative difference is still only 
0.02. We note that observed C epheid rotation rates (y rot sin i) 
are usually less than 20 km s" 1 dNardetto et al. II2006I) . 

4.7. Numerical Resolution 

Finally, we have tested the effects of the numerical resolu- 
tion of the models. In this work we have used 120 zones in 
our models, which is sufficient to generate decent light curves. 



As the phase lag itself is sensitive to the structure of the outer 
stellar layers, different zoning in this region can potentially 
influence the lag between the velocity and luminosity max- 
ima. We have recomputed some of our Galactic models with 
80 and 300 zones. The linear phase lags turned out to be rather 
insensitive to the zoning. The difference is the smallest for the 
low period models, and increases toward higher masses. The 
maximum deviation is ^5 0.02. Therefore we conclude that 
different model zoning has practically no influence on A$i . 

5. PULSATION AMPLITUDES 

So far we have compared bolometric theoretical light curves 
to observed V magnitudes when constructing the phase lags. 
In principle, this inconsistency could cause a problem. Obvi- 
ously the velocity curves are free from this caveat. However, 
we think that V phase lags are close to their bolometric coun- 
terpart, because the T e g of the Cepheid variables ensures that 
the maximum of their spectral energy distribution is close to 
the V band. The values f or V and I, for exa mple, are 
indeed found to be very close dNgeow et al. Il200"3l) . For a 
more detailed comparison of the V, B, R, I and bolomet ric 
light curve parameters we refer to Simon AMoffetiOHl. 

However, in order to compare our computed amplitudes to 
the observed Galactic amplitudes we need to convert from our 
bolometric amplitudes to V a mplitudes: My = Mj, \ +BC. For 
this transformation we follow Kovacs (2000), viz. 

BC = 2.0727Ar-8.0634(Ar) 2 (2) 
Ar = Logr eff -3.7720 (3) 

These fits to the stat ic atmosphere models were derived for 
Beat Cepheids, but (Bea ulieu et al. 1 1200 lh obtained similar 
formulae for a broader mass range, the difference being only a 
small systematic shift between the two (ALogL = 0.02). Ne- 
glecting the logg term in the BC formula doesn't affect our 
conclusions. The introduced error is 0™1 for our lowest and 
highest temperature models, generally being much less than 
that. 

On Fig. [7] the A y Fourier coefficients of the V light curves 
are shown for the Ogloz aTf al. I ([2000) sample together with 
our full amplitude mode ls. The lowest amplitude O l Cepheid 
on the figure is a UMi dMoskalik & Og loza 2000), included 
here to complement the observed sample. F Cepheids are de- 
noted by squares, overtones by triangles. The A\ amplitudes 
are plotted for our model sequences with 4M Q to 9.5M Q for 
F, and with 4 to 7M for Ol models, computed by using the 
Resonance M-L relation and set A of convective parameters. 
Since we have not computed nonlinear model pulsations at 
very low amplitude, we have extrapolated our results to zero 
amplitude, where applicable, i.e., at the F red edge and at the 
Ol blue edge. 

Fig. [7] shows good agreement between the model calcula- 
tions and the observations. The only small discrepancy is that 
the Ol Cepheids do not seem to extend to quite high enough 
periods, but this partially is due to the coarse grid we used 
(100K). Because of this, for the 7.5M sequence we do not 
find stable nonlinear Ol pulsation, while by using a finer grid 
we could find a narrow temperature range of overtone pulsa- 
tion. 

It's important to notice that Fig. [7] also appears to show a 
large observational bias toward large amplitude Cepheids, es- 
pecially among the fundamental mode pulsators. This fact 
has important consequences for the observed phase lag distri- 
bution. This issue is discussed in Sec. [6] 
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FIG. 7. — Observed and theoretical values of the A[(V) Fourier amplitude, 
obtained with the Resonance M-L relation and the a parameter set A. The 
theoretical F models are for our sequences with masses from AMq to 9.5Mq. 
The overtone sequences which are located at the left lower side, correspond to 
sequenc es with masses rangin g from 4Mq to 7Mq . Observed A i (V ) values 
from the Oeloza et al. ( 2000) sample are overlaid for the F (squares) and Ol 
(triangles) Cepheids. 

We now turn to the effects of resonances on our amplitude 
diagram. The action of the P2 : Pq = 1 : 2 resonance between 
the F mode and the second overtone, which is at the origin 
of the well known Hertzsprung progression, is quite visible 
near Pq = 10 in the M=6M Q sequence. It causes a dip in the 
A\ values. Our models duly follow the observed amplitude 
minimum. 

A corresponding 'overtone' resonance between the fourth 
and first overtone P4 : P\ = \ : 2 appears at Pi = 3?5 — 4?0 in 
the Ol sequences of the M=5.5M Q sequence. In this case one 
has to keep in mind that it is easier to locate the resonance 
position from the velocit y curves than form the light curves. 
iFeuchtinger et al. 1 (|2000) found the resonance to be centered 
on 4?2 (although with a slightly different M-L relation). The 
small distortion and the larger dip between the M=5.5 and 
M=6M overtone amplitude curves corresponding to this res- 
onance are not seen in the observed sample. The likely reason 
is the small number of observed Ol Cepheids, and the uncer- 
tain mode identification in Galact ic Cepheids. 

The Pi, : Pp = 1 : 3 reson ance dMoskalik & BuchlerlfT98l 
lAntonello & MoreThl 1 19961) has no visible effect, neither 
on the computed phase l ag nor on the A\(V) diagrams. 
Moskalik & Buchlen(ll990|) showed that anothe r possible res- 
onance, Pi : Pq = 2 : 3 (lAntonello & Morelhl 1 1996b has no 
effect unless it destabilizes the F limit cycle, in which case 
it causes the appearance of slightly alternating pulsation cy- 
cles (period doubling). Such pulsations were reported by 
Mosk alik & Buchierl (1 199 lb in radiative Cepheid models, but 
not in the sequences of turbulent models that we have com- 
puted here. No observational evidence has yet shown up to 
corroborate the existence of alternating cycles in Cepheids of 
this period range. 

To conclude, we have shown that the amplitudes of our 
models reproduce the observed amplitudes, as well as the res- 



onance characteristics. Furthermore, we have pointed out a 
deficiency of observed small amplitude Cepheids on our am- 
plitude diagram Fig. [7] These features are important regard- 
ing the following section, in which the observed and theo- 
retical phase lags are compared, and an explanation is pro- 
posed for the observed phase lag distribution for both F and 
Ol Cepheids. 

6. COMPARISON WITH OBSERVATION 
6.1. The fundamental mode phase lag 

At this time only phase lags fro m Galactic Cepheids 
have been reported in th e literature dOgloza et al. 1 120001: 
iMoskalik & Oglozall2000l) . In Fig. [5^ we have plotted these 
observational Galactic phase lags. The F Cepheids and Ol 
Cepheids denoted by squares and triangles, respectively. They 
are superposed on the linear and nonlinear phase lags that 
we have computed with our 'standard' reference parameters 
in this paper, namely the convective parameter set A ( a.) 
panel), the Resonance M-L relation and a Galactic metallic- 
ity (Z = 0.020). 

In Figs. [3^ and [5] it is seen that our model sequences ex- 
tend over a much broader range of periods than the obser- 
vational data. The simple reason is that our sequences are 
not constrained by the evolutionary tracks (except for the M- 
L relation), and thus they cover a much larger range of stellar 
parameter values (L, M, T e g). 

In addition, the observed stars strikingly occupy only the 
upper section of the region in the A$i vs. Log P diagram that 
is permitted by the modeling. This region corresponds to the 
inside region of the IS; in other words the observed F sample 
seems to avoid the vicinity of the F red edge (Fig. [3^). We 
argue that this arises from an observational amplitude bias. 
Fig. [7] suggests that the low amplitude Cepheids are missing 
from the observed sample, most likely because of observa- 
tional selection effects. We note though that there is also a 
theoretical reason for the a bsence of low amplitude Cepheids 
(Buchler & Kol llthll2002h . In all our Log P - A$, figures, 
Figs.[3t and|4]the bottoms of the downward legs of the funda- 
mental A$i curves correspond to the FRE and thus to mod- 
els with pulsation amplitudes that vanish as the FRE is ap- 
proached. In addition, towards the red edge the phase lags of 
the F Cepheid models are extremely sensitive to temperature 
as witnessed from the spacing of the dots in Fig. [3^. 

Fig- ju xtaposes the the oretical nonlinear F phase lags 
and the Ogloza et al. I (120001) sample. As discussed in §4 the 
F models sequences terminate at finite amplitude at the non- 
linear fundamental blue edge (NFBE), because the model se- 
quence switches either to double mode or to Ol single mode 
pulsation at that point. 

Overall the agreement of the model phase lags with the 
observational ones is good (differences of 0.1 rad translate 
to only about 6°). However, a few systematic discrepancies 
stand out. For example, the observational F Cepheid phase 
lags fall slightly above the theoretical ones in the middle range 
of periods. But, from a comparison of Fig. [5^ and Fig. [6] we 
see that this discrepancy would be mitigated or even removed 
if we decreased the metal content that we have assumed for 
the Galactic Cepheid model s, in line with a recen t revision of 
the solar metal abundances ( Asplu nd et al. Il2005h . 

In Fig. [HJ) lines of pulsation equi-amplitude (cutoffs) are 
overplotted from 0™0 at the FRE to 0™7. To avoid confusion 
we caution that these lines refer only to the rising parts of 
the A$i curves near the FRE. (There may be multiple equi- 
amplitude lines because the pulsation amplitudes are not just 
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FIG. 8. — Nonlinear (thick lines) phase lags A<3?i, obtained with he Resonance M-L relation, the a parameter set A, and a Galactic compositio n (bottom mass 
labels), a.): first overtone Cepheids, b.): fundamental Cepheids. The observational phase lags of the Galactic Cepheids I Ogloza et al. 2000) are denoted by 
triangles and squares, respectively. Individual Ol stars are discussed in the text. On the left panel, dashed lines shows models computed with the Girardi M-L 
relation (top mass labels). On the right, V magnitude cutoff lines (dashed lines) help to define the region of normal amplitude Cepheids (from 0™0 to 0™7). The 
FRE and the OBE are shown for reference. 

Girardi M-L relations dip strongly below the Resonance M-L 
relations. 

In Fig. [8^ the Ol observational points seem to fall a little 
higher than the ones computed for Galactic models. Again, 
reducing the metallicity would at least partially mitigate the 
difference. Most of the discrepancy occurs because of 3 re- 
ally egregious stars, two with $i > -0.3 (IR Cep and V950 
Sco), and the long period V440 Per with $i < -0.7. For 
the first two, strong evidence has been given on the basis of 
the Fourier decomposition parameters of the light and radial 
velocity curves that they are Ol pulsators (Moskalik, private 
communication). However, according to their location in the 
phase lag diagram, Fig.[3k, at least, they would quite naturally 
fit in the F, rather than in the Ol group. Alternatively, the 
disagreement for IR Cep would also disappear if its metal- 
licity was somewhat below the average Galactic metallicity, 
although this would not apply as well for the longer period 
V950 Sco. 

V440 Per has a low am plitude of Ai=0.09 5, and an enor- 
mous error bar for A$i (Ogl oza et al. | |2000). It is classified 
as the longest period Galactic Ol Cepheid. One notes that the 
Ol Cepheids MY Pup and its immediate neighbors to the left 
fall right in the middle of the resonance region. The employed 
M-L relation has an effect on the location of this resonance, 
and this in turn affects the results for A$i. A comparison 
between the phase lags in Fig. [8^ with the Girardi and the 
Resonance M-L relations suggests that an M-L relation that 
lies a little higher than the Resonance M-L in Fig. [T]in the 
5.5 to 6.0 Mq range should push the computed A$i a little 
further to the right and improve the agreement for the stars 
with Logfi > 0.6, and in particular for MY Pup. However, it 
is unlikely that this would solve the problem with V440 Per, 
but we note that the phase lag of this star is compatible with 



a function of period and A$i; the same amplitude can occur 
more than once along the same arch, but generally not on the 
same line.) 

Based on Fig. HJ) we can see that the lower limit of the 
observed amplitude distribution is well represented by the 
m 7 and m 3 cutoff lines for P£ lO^and P<> I0 d , respectively. 
This dichotomy and the sharp change in the equi-amplitude 
lines on the A$i - log P plane around 10 days has its origin 
in the P2 : Po = 1 : 2 resonance discussed above. Fig. to- 
gether with Fig. [7] confirm our hypothesis that observational 
selection effects make the observed sample gravitate toward 
the high amplitude region. 

6.2. The first overtone phase lag 

On Fig. [8^ we display the region of Ol models separately, 
to avoid overcrowding. Five individual stars that are men- 
tioned below are labeled and plotted with open symbols. As 
discussed in §4 the Ol models sequences terminate at finite 
amplitude at the NORE, because the model sequences con- 
tinue as F single mode pulsation to the right of that point. The 
P4 : Pi = 1 : 2 resonance is seen to have a remarkable effect 
on the M=5.5M and M=6.0M Q phase lags. As is seen most 
sharply in the 5.5M Q sequence the resonance causes A$i to 
drop sharply below its value at the blue edge as the nonlinear 
overtone red edge (NORE) is reached. 

In Fig. [8^ we have also shown the phase lags of 4.5, 5.0, 5.5, 
6.0 and 6.5M sequences computed with the Girardi M-L re- 
lations. They are shown as dashed lines. (The 6M Q curves are 
almost indistinguishable for the two different M-L relations). 
Interestingly we find that while the low mass curve is shifted 
to considerably smaller periods, the envelope of the accessible 
A$i curves is affected only slightly. (The reason for the shift 
to lower periods is of course consistent with Fig. [T] where the 
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the theoretical ones for small amplitude F Cepheids as seen 
in Fig. |5h, suggesting that perhaps it is a small amplitude F 
Cepheid rather than an overtone Cepheid. 

a UMi is a well-known Cepheid, most probably pulsating 
in the first overtone mode. This star has been showing dra- 
matic amplitude changes in the last decades dTurner et al. I 
120051) . Currently being a l ow amplitude pulsator with Ay = 
0! n 015 dFernie et al. 1119951) . it resides re latively close to our 
blue edge on the Log P - A$; plot (Moskalik & Ogloza 
2000). However this agreement should be taken with a grain 
of salt, because this star can be in the phase of cros sing the in- 
stability strip for the first time dTurner et al. f 2005) thus bear- 
ing a different M-L relation tha n most of the Cepheids (see 
e.g., Fig. 2 of iBuchler & Szabol (120071 )). If this is the case, 
then evaluating it's position on Fig. [8^ would require more 
computational work with the appropriate M-L relation, which 
is out of scope of this paper. 

Finally we note that on a large s cale the Ol phase la g distri- 
bution appears essentially linear (Oglo zaef al. 1 12000'). If the 
three egregious stars are removed from the Ol Cepheids then 
not only does the observational phase lag distribution appears 
curved, but that our results closely match its general shape 
(albeit being a little lower). We have also shown that the cur- 
vature of the distribution has its origin in the Pi : Pq = 2:3 
resonance. 

6.3. The phase lag of other radially pulsating Cepheids 

Recently ultra-low amp litude (ULA) Ceph eids have been 
discovered in the LMC ( IBuchler et al. 1 120051) . Based on 
this work we predict that ULA Cepheids pulsating in the F 
mode and close to the fundamental red edge, would feature 
a greater (negative) phase lag, and would fall below the cur- 
rent observed range of the regular amplitude Cepheids. Also, 
ULA Ol Cepheids with periods far from the resonance center 
would be close to the lower envelope of the predicted phase 
lag distribution. 

Our Fig. [5] offers a method for radial mode identification 
for Cepheids that is additional and complementary to the 
usual relative Fourier decomposition coefficients (/?2i,^2i> 
etc). We see that it offers not only a criterion for the F and 
Ol Cepheids, but also for the elusive Galactic 02 Cepheids. 
There clearly is an observational bias against low period - low 
amplitude Cepheids. However, considering the existence of a 
Galactic 01&02 beat Cepheid (CO Aur) it is quite likely that 
pure 02 Cepheids exist as well. The phase lag has the po- 
tential of sorting the elusive Galactic 02 Cepheids from Ol 
Cepheids. 

Observed and theoretical Fourier parameters of Cepheid 
light curves have be en often investigated, see e.g., 
Simon & Kanburl d 19951) for lon g perio d F Cepheid Fourier 
decomposition, ISimon & Davisl ([1983) for F phase lags. 
Other previous theoretical work s on Ol Cepheid pulsa- 
tions dAntonello & Aikawa|| 19951) and on 02 Cepheid pulsa- 
tions (jAntonello & Kanb url 19971) Cepheids also discussed the 
phase lags. However, these authors ignored convection and 
produced purely radiative models, an approximation which 
makes sense only near the blu e edge. The resulting large 
range of F and Ol phase lags led Antonello & Aikawa ( 1995) 
to conclude that A$i is not a useful discriminator for the 
pulsational mode. One problem may be that radiative mod- 
els have excessively large amplitudes. Because the nonlinear 
phase lags depend largely on the amplitude, as we pointed out 
in Sec. 14.41 we think that to address the phase lag problem 



correctly, convection can not be left out from the models. 
7. CONCLUSION 
In this paper we have reexamined the classical Cepheid 
phase lags, A$i = - <fi™" 8 , with our hydrocode which in- 
cludes a time-dependent mixing length model for turbulent 
convection. 

As expected, the linear and nonlinear (full amplitude) phase 
lags approach each other when the pulsation amplitude is low, 
as for example near the fundamental red edge of the instabil- 
ity strip. The nonlinear phase lags generally fall below the lin- 
ear ones except for high period models where they lie above. 
The differences between the linear and nonlinear phase lags 
(| A(A$i)|) are generally at most 0.20 for our Galactic mod- 
els with the best convective model parameters (set A). In case 
of lower metallicities, the |A(A$i)| can be slightly larger. 

In the linear regime we find that high mass models of 
low metallicity can have positive A$j. However, for Galac- 
tic composition the nonlinear A$i are always found to be 
negative, i.e., the radial velocity always lags the luminosity 
slightly. 

There is a small metallicity effect: Lower metallicity mod- 
els pulsating in the fundamental mode show phase lags that 
are slightly higher than the Galactic values. There are how- 
ever no sufficiently good observations to confirm this. 

The convective a parameters have a strong influence on the 
modal selection in the instability strip, as well as on the pul- 
sation amplitude, which is then reflected in the Period - A$i 
diagram. Similarly, the M-L relations have only an indirect 
effect on the phase lags through a modification of the mass- 
period relation. 

We obtain a good general agree ment between the obse rved 
Galactic Cepheid phase lag data of [Qgloza~e f al. (200(3) and 
our models. For the F Cepheids this agreement can be fur- 
ther improved when the observed amplitude distribution is 
taken into account. By applying cutoff magnitude (equi- 
amplitude) lines we can restrict the model phase lags region 
to the approximate observed range [-0.2; -0.4]. An almost 
perfect agreement can be achieved by lowering the assumed 
metallicity of Z = 0.0 20 in accordance with the suggestion of 
lAsplund et~aT\ (120051) . 

Our computed Ol phase lags seem to fall on a steeper re- 
lation than the observations, but this discrepancy would be 
mitigated or removed if we could reclassify two stars, IR Cep 
and V950 Sco as F Cepheids (even though Fourier decom- 
position parameters of the light and radial velocity curves in- 
dicate otherwise). Furthermore, the phase lag of the outlier, 
V440 Per, which is normally classified as an overtone, would 
be compatible with the theoretical phase lag of a low ampli- 
tude F Cepheid. With the omission of these three stars our 
results match the shape of the observational distribution of 
phase lags. 

We have shown that resonances have an impact on the phase 
lag. The P% : Pq = 1 : 2 resonance which is present in the fun- 
damental mode Cepheids at around lO'' period and causes the 
well-known Hertzsprung-progression, has a small effect on 
the minimum in the nonlinear phase lags (Fig. [3^ and Fig.|5^). 
As for the Ol models, the P4 : Pi = 1 : 2 resonance generates 
an abrupt widening of the 'either F or Ol' region at around 
the resonance center (Szabo et al. 2007, in preparation). In 
this region A$i is more negative than for models computed 
far from the resonance, and this effect seems to have great 
influence on the observed O l phase lag dis t ributio n. 

We confirm the finding of Ogloza ef al. I d2000t) and lSimonl 
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(1984) that the A$i is a good indicator of the pulsational 
mode that is independent of and supplemental to the usually 
used relative Fourier decomposition coefficients (/?2i,^2i> 
etc). In principle, it is also a good criterion for identifying sec- 
ond overtone pulsators. However, no 02 Cepheids have been 
detected in the Galaxy, but there are a number of candidates 
in the Magellanic Clouds. For an independent confirmation 
of their status as 02 on the basis of the phase lag it would be 
necessary to acquire radial velocity data on these objects. 

Ultra-low amplitude Cepheids should be easily discernible 
based on their phase lags. Our work predicts that ULA 
Cepheids pulsating in the fundamental mode and close to the 
fundamental red edge, would feature a greater phase differ- 
ence, and would lie below the observed range [-0.2; -0.4] of 
the large-amplitude Cepheids. ULA first overtone Cepheids 
(at the blue edge) would be similarly separable, except in the 
vicinity of the 'overtone' resonance. Therefore obtaining ra- 
dial velocity data and existing light curve parameters could 
further prove their nature, i.e., radial pulsation. Unfortunately, 
getting V r curve for a small amplitude Cepheid in the MCs re- 



mains a huge technological challenge. 

Our survey has produced concurrently a rather complete 
theoretical overview of the nonlinear Cepheid instability strip, 
which delineates the regions of F, 1 , 02 pulsation, of F&O 1 
and 01&02 beat pulsation, and regimes of hysteresis, such 
as the 'either F or Ol' regions and 'either Ol or 02' re- 
gions. These findings w ill be presented in a separate paper 
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